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The electronic band structure of CeAgSb2 has been calculated using the self-consistent full poten- 
tial nonorthogonal local orbital minimum basis scheme based on the density functional theory. We 
investigated the electronic structure with the spin-orbit interaction and on-site Coulomb potential 
for the Ce-derived 4f orbitals to obtain the correctly ground state of CeAgSb2 . 

PACS numbers: 71.10.Hf, 71.18.+y, 72.20.Eh, 75.30.Mb 



I. INTRODUCTION 

The tenary compounds CeXSb2 (where X=Au, Ag, 
Cu, Ni, or Pd) have attracted considerable attention 
recently because they exhibit a large variety of physi- 
cal properties such as heavy fermion, Kondo insulating, 
anisotropic transport and magnetic ordering behavior. 
The various physical properties depend on the hybridiza- 
tion between f electrons and conduction electrons, which 
is characterized by a Kondo temperature Tk- With a 
small hybridization the system exhibits a local magnetic 
moment and orders magnetically, which can be described 
by the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action. By increasing the hybridization, the Kondo effect 
increases and the ordered magnetic moment decreases 
gradually. With further increase of the hybridization, 
the system will go into a heavy fermion or intermediate 
valence regime. 

CeAgSb2 is a Ce-based Kondo lattice system with low- 
temperature ferromagentic ordering. Sologub et al. [l| 
reported this compound has weak ferromagnetic oder in 
polycrystalline samples with a net ferromagnetic moment 
of 0.15 fis/Ge at 5 K. Several different groups investi- 
gated the magnetic properties of CeAgSb2 with conflic- 
ing results for the magnetic ground state. HHQ. In 
particular, for single crystal smaples, the magentization 
of CeAgSb2 is anisotropic 0. In this case the magneti- 
zation for the magnetic field applied parallel to the c- 
axis shows a ferromagnetic order with a moment of 0.37 
fis/Ce above 0.025 T at 2 K. For H±c, the magneti- 
zation increases nearly linearly to 1.2 /x^/Ce below 30 
kOe and then remains nearly constant for higher fields, 
possibly indicating the presence of a metamagnetic tran- 
sition. Above 200K, the inverse susceptibility is linear to 
the Curie- Weiss law giving an effective moment of 2.26 
Hb I Ce. In zero field, the resistivity of CeAgSb2 increases 
rapidly from 1.16 /ificm at 2K to a maximum of 88.1 
/ificm at 18.2 K. At 9.7K, a sharp change in the slope 
of the zero field resistivity is observed, consistent with a 
loss of spin-disorder scattering associated with the mag- 
netic ordering as well as the possible suppression of the 
Ce hybridization due to the ferromagnetic ordering. At 
higher temperatures, the resistivity determines a broad 
local minimum near 150K. The temperature-dependent 
resistivity of CeAgSb2 is typical of a Kondo lattice sys- 
tem. Preliminary inelastic neutron scattering measure- 



ments indicate a Kondo temperatur Tk between 60 and 
80K Q. The magnetoresistence of CeAgSb2 is either 
positive or negative, depending on the temperature and 
orientation of the applied field. This complexity is due 
to competing contributions from the Ce hydridization, 
magnetic ordering, and the electronic structure of the 
compound. The magnetoresistence of a paramagnet is 
expected to be negative, since the local moments try to 
align parallel to the applied field, which essentially re- 
duces the spin disorder scattering as the field is increased. 

In order to understand the electronic and magentic 
properties of CeAgSb2, we need the electronic band 
structure studies based on the density functional theory. 
In this work, the precise self-consistent full potential local 
orbital minimum basis band structure scheme (FPLO) 
are employed to investigate the electronic and magnetic 
properties of CeAgSb 2 with LDA, LDA+U and fully rel- 
ativistic shemes. We consider the effect of magnetism on 
the band structure and compare with experiment. 



II. CRYSTAL STRUCTURE 

CeAgSb2 crystallizes in the primitive tetragonal 
ZrCuSi2 structure, which consists of Sb-CeSb-Ag-CeSb- 
Sb layers along [001] direction. This crystal structure is 
described in detail by Myers et al. [f|. It belongs to the 
P4/nmm space group with Ce occupying the 2c site, Ag 
the 2a, and Sb the 2b, 2c sites. The sites are given, in 
units of (a, a, c): Ce(0.25, 0.25, 0.26), Ag(0.75, 0.25,0.5), 
Sb(0. 25,0. 25,0.84), Sb(0.75, 0.25, 0). We used experi- 
mental lattice constants, a,— 4.36 A and c=10.41 A, in 
the calculation described below. There are two formula 
units in the primitive cell. 



III. METHOD OF CALCULATIONS 

We have applied the full-potential nonorthogonal local- 
orbital minimum-basis (FPLO) scheme within the local 
density approximation (LDA). 0] In these scalar relativis- 
tic calculations we used the exchange and correlation po- 
tential of Perdew and Wang.0 Ce 4s, 4p, 3d, Ag 5s, 5p, Ad 
and Sb 4s, Ap states were included as valence states. All 
lower states were treated as core states. We included the 
relatively extended semicore 4s, 4p, 4d, 4f, 5s, 5p states 
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FIG. 1: The full LDA nonmagnetic band structure of 
CeAgSt)2 along the symmetry lines. 

of Ce and 4s, 4p states of Sn and Ag as band states be- 
cause of the considerable overlap of these states on near- 
est neighbors. This overlap would be otherwise neglected 
in our FPLO scheme. Ce 6p states were added to in- 
crease the quality of the basis set. The spatial extension 
of the basis orbitals, controlled by a confining potential 
(r/ro) 4 , was optimized to minimize the total energy. The 
self-consistent potentials were carried out on a k mesh of 
18 k points in each direction of the Brillouin zone, which 
corresponds to 550 k points in the irreducible zone. A 
careful sampling of the Brillouin zone is necessitated by 
the fine structures in the density of states near Fermi 
level E F . 



IV. RESULTS AND DISCUSSION 

We first show the full band structures of CeAgSb2 
within LDA scheme in Fig. 2] The Sb 5s bands lie be- 
tween -11.5eV and -6.5eV. Between -6.5 and the Fermi 
level there are of mixed Sb 5p and Ag 4d states. Those 
very flat bands near the Fermi level are mainly the Ce- 
centered 4f characters. A prominent feature of the band 
structure near Ep, besides the 4f bands, is the Ce 5d 
character which hybridizes with the Ce 4f bands. We 
also study the on-site atomic-like correlation effects be- 
yond LDA by using LDA+U approach in a rotation- 
ally invariant, full potential implementation 8]. Minimiz- 
ing the LDA+U total energy functional with spin-orbit 
coupling(SOC) treated self-consistently generates not 
only the ground state energy and spin densities, but also 
effective one-electron states and energies that provides 
the orbital contribution to the moment and Fermi sur- 
faces. The basic difference of LDA+U calculations from 
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FIG. 2: Top panel:The LDA band structure of CeAgSb 2 
along symmetry lines. The very flat bands near the Fermi 
level are the Ce 4f bands. Middle panekThe band structure 
within the LDA+U scheme showing that the 4f bands are split 
in a 1-3-3 fashion from bottom up. Bottom panekThe fully 
relativistic band structure of CeAgSb2 showing that spin-orbit 
coupling splits the 4f states into two manifolds. 
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FIG. 3: Projected density of states of CeAgSb2- Top panel: 
total, plotted positively for majority, negatively for minority. 
Bottom panel: projection of the Ce 4f and Sb 5p, showing that 
Ce 4f character dominates the states near the Fermi level. 



the LDA is its explicit dependence on the on-site spin and 
orbitally resolved occupation matrices. The Coulomb 
potential U and the exchange coupling J for the Ce 4f 
orbitals have been chosen to be 7.5 and 0.68 eV respec- 
tively. The resulting band structure calculated within 
the LDA+U scheme is shown in middle panel of Fig. 
[21 We observe that the crystal field splittings of Ce 4f 
bands within LDA are quite small and in fact difficult to 
identify due to hybridization with itinerant bands. From 
LDA+U, the 4f bands are still very flat but are split (in 
a 1-3-3 fashion from bottom up) by some combination 
of the crystal field and the anisotropy of the U interac- 
tion by a total 1.6 eV. We also calculated fully relativis- 
tic band structure to see the spin-orbit coupling effects, 
which is shown in the lower panle of Fig. [21 As ex- 
pected spin-orbit coupling splits the 4f states into two 
manifolds, located 0.4 eV and 0.05 eV above the Fermi 
level, the 4f 7 / 2 and the 4f 5 / 2 multiplet respectively. 

The atom and symmetry projected densities of states 
(PDOS) shown in Fig[3] clarify the characters of the 
bands. Because only the DOS distribution near the Fermi 
level determines the magnetic properties, we concentrate 
our attention upon the DOS in the vicinity of the Fermi 
level. At this range, the valence states for Ce or Sb 
atoms are dominated by 4f and 5p electrons, respectively, 
and the contributions from other electrons are negligibly 
small. Because of hybridization, they distribute in a wide 



energy range and extend to the unoccupied states, above 
the Fermi level. These states remain almost unchanged 
when the spin-orbital interaction or on-site correlation 
potential are taken into consideration. So only the DOS 
obtained by LDA calculation for Ce-4f and Sb-5p have 
been plotted. From the Fig. [3J it can be found that 
within the LDA calculation, the DOS at the Fermi level 
are larger and they are mainly of Ce-derived 4f states. 
The contribution from Ag-derived 4d or Sb-derived 5p 
states is smaller. When the spin-orbital couplilng along 
the axis (0 1) is taken into account, the 4f orbitals are 
slightly modified. The spin 4f states become wider and 
the energy shift between centers becomes larger. This 
is due to the partial splitting between the degenerate 4f 
states. 

Band calculation! with the LDA framework cannot 
yield correctly magnetic moment for many Ce compounds 
because of the strong correlation interaction between f 
orbitals. The spin-orbit interaction in these systems is 
sometimes large and the orbital contribution to magnetic 
moment cannot be neglected. The calculated magnetic 
moment of CeAgSb2 within LDA scheme is 0.61 /i^/Ce. 
it is mainly form the Ce-derived 4f orbitals, while almost 
no contribution from Ag. This is in agreement with the 
fact that the transition metal T, except Mn, does not 
carry magnetic moment in CeTSb2 compounds. When 
the on-site correlation potential is added to the Ce 4f 
electron, the degeneracy between the different f orbits 
would been lifted and the Hund's rules dominate the lo- 
cally occupied 4f electrons, which yields the total mag- 
netic moment 0.54 /j^/Ce. With fully relativistic scheme 
we calculated the magnetic moment of the value 0.55 
/is/Ce. As mentioned above, several groups reported the 
experimental studies on the magnetic moments which do 
not agree with another, so further experimental investi- 
gations are required. 

Density functional calculations are very reliable in cal- 
culating the instability to ferromagnetism. The pres- 
ence of an electronic instability is signaled by a diver- 
gence of the corresponding susceptibility. In the follow- 
ing we study the uniform magnetic susceptibility using 
the method of Janak ^(j ■ The uniform magnetic suscep- 
tibility of a metal can be written as 



A l-N(E F )r y ' 

where the numerator stands for the Pauli susceptibil- 
ity of a gas of non-interacting electrons proportional to 
the density of states at the Fermi level N(E F ), and the 
denominator represents the enhancement due to electron- 
electron interaction. Within the Kohn-Sham formalism 
of density functional theory the Stoner parameter / is re- 
lated to the second derivative of the exchange-correlation 
fuctional with respect to the magnetization density. We 
have evaluated, within the density functional theory for- 
malism, the Stoner enhancement of the susceptibility 
X = i-in\e f ) = S X°> where = 2fi 2 B N(E F ) is the 
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non-interacting susceptibility and S gives the electron- 
electron enhancement in terms of the Stoner constant I. 
We have calculated I using both the Janak-Vosko-Perdew 
theory [l^ and fixed spin moment calculations . The 
calculated density of states and Stoner parameter nor- 
malized per are N(Ep)=5.8 states/eV and 7=0.19 eV. 
This gives IN(Ef) = 1-1, larger than unity, correspond- 
ing to a ferromagnetic instability. 

Heavy fermion compound is characterized by a larger 
electronic specific heat coefficient 7. CeAgSb2 is a moder- 
ate heavy fermion compound with 7 = 65 mJ/K 2 mol|l2j. 
The large specific heat coefficient of CeAgSb2 compound 
could not be yielded by our band calculation. This can 
be seen from the calculated electronic structure. It can 
be found that the toal number of DOS at the Fermi 
level is about 5.8 states/eV, which corresponds 7;, = 13.6 
mJ/K 2 mol and underestimate the experiment value by 
a factor of 4.8. The discrepancy between the band cal- 
culation and experiment for specific heat coefficient is 
attributed to the formation of quasiparticle. There is ex- 
change interaction J between the local f and the conduc- 
tion electrons in CeAgSb2. The ground state of Ce com- 
pound is determined by the competition of the Kondo 
and indirect RKKY interaction. With a large J, the 
Kondo coupling becomes strong and the system located 
at the borderline of magnetic-nonmagnetic transition. 
The exchanging interaction between the local f electron 
and the conduction electrons will result in the forma- 
tion of quasiparticle. It has a larger mass compared with 
bare electron and the enhancement of mass increases with 
the increase of exchanging. Because of the volume con- 
traction, the exchange interaction between the f and the 
conduction electrons is large in CeAgSb2. This will re- 



sult in the f electrons to behave like itinerant electrons 
and the narrow f bands to be located at the Fermi level. 
On the other hand, when the exchanging interaction be- 
tween f and conduction electrons is smaller, the occupied 
4f orbitals are located near the Fermi level while the un- 
occupied 4f orbitals are at the conduction bands. The 
quasiparticle mass is appropriate to the number of DOS 
at the Fermi level. So the quasiparticle mass is largely 
enhanced in CeAgSb 2 . Indeed, it has been shown that 
when the Ce 4f electrons in CeAgSb2 are treated as local- 
ized electrons, the quasiparticle mass was only enhanced 
over the band calculation by a factor 4.8. 



V. SUMMARY 

In this article we showed the results of three different 
electronic band structure calculations. It shows that the 
Coulomb potential on Ce 4f orbitals and spin-orbit inter- 
action is a key factor to understand the electronic and 
magnetic properties of CeAgSb 2 . When the Coulomb 
potential is added to the Ce 4f orbitals, the degeneracy 
between the different f orbits would be lifted and they 
are split into lower Hubbard bands at the Fermi level 
and unoccupied upper Hubbard bands in the conduction 
band. The exchange interaction between local f electrons 
and conduction electrons play an important role in the 
heavy fermion characters of them. And the fully relativis- 
tic band structure sheme shows that spin-orbit coupling 
splits the 4f states into two manifolds, the 4f//2 and the 
4f5/2 multiplet. 
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